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Abstract. We study the dynamics of Bose-Einstein condensates in accelerated 
c/3 two-dimensional optical square lattices by numerically solving the Gross-Pitaevskii 

equation. We consider the regime with negligible mean-field interactions and examine 
in detail the pulses of atom clouds ejected from the condensate due to Landau-Zener 
tunnelling. The pulses exhibit patterned structures that can be understood from the 
momentum-space dynamics of the condensate. Aside from conceiving realization of 
a pulsed two-dimensional atom laser, we demonstrate that, by exploring the band 
structure of the lattice, Landau-Zener tunnelling and Bragg reflection of the condensate 
C$ inside the optical lattice can provide means for manipulation of coherent atom waves. 
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1. Introduction 

The past two decades have witnessed the fruitful interplay between condensed matter 
physics and ultracold atom physics. While cold atom systems have provided new arena 
for condensed matter physics, at the same time, concepts and ideas from condensed 
matter physics have opened up new directions for cold atom physics. For instance, by 
loading ultracold atoms into optical potentials due to standing laser waves (i.e. optical 
lattices), it is possible not only to build up "quantum simulators" for condensed matter 
systems (see, for example, jl]), but also to control and manipulate cold atoms utilizing 
concepts from traditional solid-state systems (see, for instance, [2]). In this work, we 
wish to further explore the latter lines of investigation. 

Electrons in solids are subject to periodic potentials due to the array of ion cores 
that constitutes the lattice. In ideal situations, the electron eigenfunctions take the 
form of a particular structure that is called the Bloch wave, and the corresponding 
eigen-energies cluster into energy bands that are separated by energy gaps (see e.g. [3]). 
As a consequence, electrons in solids can behave quite differently compared with those in 
free space. For instance, when a weak uniform static electric field is applied to the solid, 
instead of being uniformly accelerated, the electrons would carry out oscillatory motions 
(the Bloch oscillation; see e.g. jl]). This oscillatory motion can be understood as Bragg 
reflections of the Bloch states in momentum space (or fe-space) at regions where energy 
gaps occur. If the bias field is large, the electrons can gain sufficient energy that would 
lead to finite probability for tunnelling across the energy gaps, resulting in transitions to 
higher energy bands (the Landau- Zener tunnelling; see, for example, [1]). In this work, 
we wish to consider the cold atom counterpart of the above construction. Namely, we 
shall consider a cloud of Bose-Einstein condensate subject to a periodic potential due 
to standing laser waves that form an optical lattice. We shall examine dynamics of the 
condensate under the action of an external bias that corresponds to the uniform static 
electric field above. Experimentally this bias can be furnished in a number of ways: 
one could make use of direct gravity pull on the atoms, or introduce time-dependent 
frequency shifts between the lattice beams so that an effective acceleration on the atoms 
occurs (see e.g. [S1E]). For spin polarized condensates, it is also possible to generate the 
acceleration by means of magnetic field gradients [7j. 

The dynamics of ultracold atoms in one-dimensional (ID) optical lattices subject 
to uniform accelerations have been investigated extensively [2]. In addition to the 
observation of Bloch oscillations [8], Wannier-Stark ladders [9], and Landau-Zener 
tunnelling [TU] for condensates in accelerated ID optical lattices, possible realizations 
for ID "atom laser" have also been achieved [H] using such systemsjf] By loading Bose 
condensed 87 Rb atoms into a vertical ID optical lattice, Anderson and Kasevich have 
demonstrated that for appropriate lattice strengths, pulses of atom clouds can tunnel 
out of the optical lattice due to the gravity pull. The generation of pulses of coherent 

| For other realizations of atom lasers (involving and not involving accelerated optical lattices), see 
e.g. [7] and references therein. 
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atom waves thus suggests the system being an ID pulsed atom laser [HI [12] . For two- 
dimensional (2D) systems, there have also been studies of Bloch oscillations, Landau- 
Zener tunnelling (at lower energy bands; cf . below) [131 E] > an d Bloch-Zener oscillations 
[15] of ultracold atoms in accelerated 2D optical lattices|§] The richer dynamics of the 
condensates in 2D lattices has been proposed to be possible means for manipulation 
of the coherent atom waves [131 H2 US]- Along these lines of investigation, we shall 
examine in this work the Landau-Zener tunnelling of atom waves within 2D optical 
lattices. More specifically, we will be interested in physical configurations similar to 
the ID system considered by Anderson and Kasevich [TTl[T2] . To simplify the analysis, 
however, we will focus in this work on regimes where the nonlinear mean-field interaction 
is negligible (cf. [HI H2J). Effects of the mean- field interaction in related problems are 
investigated in [201 |2~T] . 

In order to conceive a 2D atom laser, we shall consider a condensate loaded into 
a uniformly accelerated 2D optical square lattice in the regime where Landau-Zener 
tunnelling is the predominant effect. In particular, since the atoms need to tunnel out 
of the lattice to form the lasing pulses, dynamics of the condensate at higher energy 
bands will be essential (cf. [131 El)- As we will see below the more complicated energy 
landscape at higher energy bands turns out to enrich the structure for the output pulses 
of the atom laser. We shall study the dynamics of the condensate by numerically solving 
the Gross-Pitaevskii (GP) equation (see e.g. [6]). For appropriate accelerations we shall 
find pulses of atom clouds tunnelling out of the condensate, which manifest a patterned 
structure. We will show that these are consequences of Bragg reflection and Landau- 
Zener tunnelling at higher energy bands. For energy states in these higher energy bands, 
unlike their ID counterparts [UJ, characteristics of Bloch waves still manifests in the 
atom waves; even though they are quite close to free-particle states. By means of k- 
space analysis, we shall relate propagation of the coherent atom waves in real space 
to the fc-space dynamics of the condensate. As we will soon recognize, aside from a 
2D atom laser with patterned pulses, our result can help envisage beam splitters for 
coherent atom waves. 



In the following, we shall start in section 2.1 with an introduction to our theoretical 



formulation. Numerical results for our simulations will then be presented in section 2.2 



In section 2J3 we will try to understand these results base on fc-space dynamics of the 
condensate. We will then demonstrate in section [3] how this analysis can help provide a 
useful tool for manipulation of coherent atom waves. A brief conclusion is supplied at 
the end of the article. 

§ There have also been related theoretical and experimental studies in photonic crystals [T6 l [T7 l fT8 l ll9j. 
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2. Two-dimensional atom laser 

2. 1 . Theoretical formulation 

Let us consider a Bose-Einstein condensate of atoms with atomic mass m in a two 
dimensional optical lattice [22J . We shall look into the regime where the low temperature 
ground-state dynamics of the condensate can be described by the GP equation (see 
e.g. ©) 

m^{ P ,t)= (-^v 2 p + v ext ( P ,t) + u\^( P ,t)A y( P ,t) , (i) 

where \& is the condensate wave function, p = (x,y) is the 2D position vector, 
Vp = (^2+^2) is the 2D Laplacian operator, V ext is the total external potential 
acting on the condensate, and U is the mean- field coupling constant. In order to examine 
dynamics of the condensate in an accelerated optical lattice, here the external potential 
Vext consists of two parts: one due to the optical lattice, and the other from a uniform 
acceleration. The optical lattice potential that we shall consider is similar to that in 



Vop t (x,y) = V [cos 2 (q x x) + cos 2 (q y y)] , (2) 

where V is a real constant, q x , q y are wave numbers for the laser beams that build up 
the optical lattice. In experiments the lattice beams usually have the same wavelength, 
so that q x = q y = q. We shall thus consider this case, so that pj corresponds to a 2D 
square lattice with lattice constant a = n/q. If the lattice is uniformly accelerated with 
acceleration g, there would be the corresponding potential 

V g (x,y) = -mg-p. (3) 

In order to maintain the phase coherence of the atom laser that we have in mind, 
it is crucial to avoid large potential strengths [12] . This is because for high lattice 
strength Vo, if one regards the condensate as a collection of individual atom clouds over 
the lattice sites, with increasing Vo the wave function overlap between neighbouring sites 
would drop and, at the same time, the mean-field interaction f/|\I/| 2 would grow due to 
the more localized wave function and stronger coupling U [23] • These effects combine 
to reduce the phase coherence among the condensates on different lattice sites [T2| [24] . 
This can be understood as a consequence of the number-squeezing of the condensate 
on each lattice site due to the mean-field interaction when Vo is large, which enhances 
the (uncorrelated) phase fluctuations among the lattice sites [24J. It is essential to 
recognize the role of the mean-field interaction in this decoherence mechanism: with 
U — 0, as long as the wave function overlap between adjacent sites remains finite, there 
would not be any decoherence even for large Vo- We will thus be interested in the 
regime where the mean- field interaction f/|\l/| 2 in (jTl) is negligible compared with the 
strength of the optical potential|j}] In experiments, this would correspond to the cases 

|j Since the mean-field interaction is proportional to the (s-wave) scattering length and the density of 
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where well-resolved Bragg peaks can be visible when the condensate is released from 
the optical lattice (Ml [25] (notice, however, the caveat pointed out by [2H]); namely 
the system is supposed to be deep on the "superfluid" side of the phase diagram. This 
would at the same time also justify our analysis of the condensate dynamics based on 
the GP equation, since for large potential strengths the system could make transitions 
into the Mott insulating phase [25]. Throughout this work, we will therefore consider 
the following GP equation 

d* / h 2 \ 

i^—= {-—V 2 p + V [cos 2 ( y qx) + cos 2 ( y qy)]-mg-p\^. (4) 

The time evolution of the condensate can then be analyzed reliably using the Crank- 
Nicholson method [29], which we shall now turn to. 



2.2. Atom laser 

Let us suppose we have initially in free space a stationary condensate which has a 
Gaussian wave function 



* M = ^ exp (-jP)' (5) 

where N is the total number of atoms, p = (x 2 + y 2 ) 1 ^ 2 , and a is the width of the wave 
packet. In order to transfer the condensate into the lattice ground state, without turning 
on V g , the Gaussian wave packet is first subject to the optical lattice potential V opt and 
set to evolve for a period of time. At time t = 0, we switch on V g [see pb] so that a static 
uniform acceleration g starts acting on the condensate. It is then our main purpose to 
examine the subsequent dynamics of the condensate. As explicated previously, to avoid 
complications from strong lattice potentials we restrict here to potential strengths of 
order of the recoil energy 

h 2 a 2 h 2 

2m &ma z 
We carry out numerical simulations for the above procedures base on the GP 
equation (El). The results for two different accelerations are presented in figure 11} where 
we have used a = 2.5a and N = 2-no 2 for the initial wave function (pi), and Vq = I.OEr 
for the optical potential fcl). The acceleration is set to have the magnitude g = | ^, 
which has been chosen somewhat arbitrarily; one only has to make sure that it would lead 
to appreciable tunnelling probability for the condensate (see later). For convenience, we 
use the following time scale for the time evolution of the condensate 

(?) _ h l[8E^ 



Tb (i) mga g]/ m ' (7) 

which is the time for a Bloch state to cross one complete Brillouin zone (namely, 
the Josephson time [11], or the Bloch period [14] ) when g is along the lattice (1,0) 

the condensate, besides tuning the lattice strength Vb (such as in [Til H21 EH ES] ) , experimentally one 
could also attain this regime by using condensates with low densities and/or reducing the scattering 
length via Feshbach resonances [5S1 [57] . 
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Figure 1. Numerical results for accelerated condensates in a 2D square lattice 
with the acceleration g along (a) the (1,0) direction and (b) the (1,1) direction. 
For both cases we have Vq = I.OEr and g = j -4, and the time scale Tb is given 
by ¥j\. Notice that in (b), for convenience, the plots are made with respect 
to the coordinates x' = (x + y)/y/2, y' = (y — x)/y/2, so that g points along 
the x'-axis. In all plots the coordinates are measured in units of the lattice 
constant a. 
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direction^] We shall henceforth adopt a coordinate system with rc-axis along the lattice 
(1,0) direction, and y-axis the (0, 1) direction. 

As shown in figure []Ja), for g along the (1,0) direction pulses of atom clouds tunnel 
out of the condensate in a manner reminiscent of that in accelerated ID optical lattices 
[ITJ [T2]. The sequence of pulses generated are equally spaced in time by T B and aligned 
entirely in the (1, 0) direction. For g along the (1, 1) direction [see figure [lib)], however, 
the pulses have a more complicated structure. There are initially three pulses generated 
from the condensate in three different directions which subsequently split into further 
sub-pulses. Clearly, this immediately brings up interesting possibilities for applications: 
besides a "patterned" atom laser, one could also use the 2D optical lattice as a beam 
splitter for coherent atom waves, as we will soon demonstrate. In the following we 
shall first examine the origin of this pattern and afterwards elaborate on its possible 
applications. 

2.3. k-space dynamics 

We shall now attempt to understand the results above from the viewpoint of fc-space 
dynamics of the condensate. This analysis will turn out very useful for real-space 
manipulations of atom waves using accelerated 2D optical lattices. 

In the absence of lattice acceleration (i.e. g = 0), the eigenstates of the GP equation 
Q for the condensate are the Bloch states (see e.g. [3]). Each Bloch state is characterized 
by the quantum numbers k (the Bloch wave- vector) and n (the band index). For the 
initial wave function (J5J), it corresponds in fe-space to a Gaussian centering at k = (0, 0) 
with a spread Sk x = 6k y = (2cr) _1 . Therefore, the subsequent time evolution of the 
condensate would yield wave packets that are superpositions of Bloch states with width 
5k ~ (2cr) -1 . For the simulations in figure n] we have a = 2.5a, thus in fe-space, instead 
of a single point, each atom cloud corresponds to a wave packet that spreads around 
a central value with area n(5k) 2 ~ -^ = ^2, which is j^ ~ 0.32% the area of one 
single Brillouin zone. In the analysis below, we shall ignore the spread in k and focus 
on the center of each wave packet, as if the atom cloud were in a sharp Bloch state 
characterized by a single Bloch wave vector. In a more quantitative treatment, one can 
take into account the spreading by averaging over the distribution of the Bloch wave 
vectors j8]. 

In view of the above, the motion of the condensate under the action of a uniform 
static acceleration g can then be dealt with using a semiclassical approach. Within each 
band (namely, for one specific n) the dynamics of the condensate follows the equations 
of motion (see e.g. [3]) 

t =i v * E - (8) 

% Notice that for g along arbitrary directions, Tg would be in general not identical to the period for 
the Bloch oscillation. In fact, it is possible to have open trajectories for Bloch oscillations [2]. In this 
case, the Bloch period does not exist. 



Manipulation of coherent atom waves using accelerated 2D optical lattices 8 

~dt =m9, (' 

where p is, as before, the 2D position vector, and E n is the eigenenergy for the Bloch 
state. As is clear from (|9]), under the action of a constant acceleration g the Bloch state 
would proceed steadily in A;-space along the direction determined by g. However, when 
the state arrives at zone boundaries where energy gaps exist, the state would have to 
either make transitions to higher energy bands (Landau-Zener tunnelling) or be reflected 
to other regions of the same zone (Bragg reflection), so that it could continue its course 
of evolution according to ^. Notice that these processes are not included in ([9]) since 
two or more energy bands can be involved. When only two states are involved, one can 
find that the tunnelling probability would be (see e.g. [3|) 

where E gap is the energy gap at the fc-point in question, and E Q is the corresponding 
free-particle energy. If the acceleration g is small, the tunnelling probability would be 
exponentially small. The state would then be Bragg reflected to another state within 
the original zone, and subsequently continue to evolve in accordance with ([8]), Q until it 
meets another zone boundary. If the state has low probability for interband transitions, 
its motion would be limited to a single zone and the motion would be oscillatory in 
both real and k spaces (which is the Bloch oscillation) [13]. However, for larger values 
of g the tunnelling probability can become appreciable. The state can then tunnel into 
higher energy bands and give rise to pulses of atom clouds in real space moving with 
group velocities determined by (J8J). As we will find out soon, the patterns of the atom 
pulses in figure [I] are consequences of Bragg reflection and Landau-Zener tunnelling of 
the atom waves. Base on this observation we will explain below, how this picture can 
be used to carry out coherent manipulations of atom waves. 

On the grounds of the foregoing analysis, we note that the state of the condensate 
migrates over the fc-space under the pull of g and hops around via Bragg reflections. 
It is thus very helpful to examine the fe-space dynamics of the condensate using the 
extended-zone picture (see e.g. [3]). Figure [2] shows the extended- zone energy landscape 
over fc-space obtained from direct diagonalization for the condensate in the 2D square 
lattice (pi). Energy gaps exist at regions where there are discontinuities in the energy 
contours, which appear as dark shading lines in figure [^b); the darker the shading is, 
the larger the energy gap is. Notice that not every zone boundary has energy gaps. 
As in the extended-zone scheme the zone index is identical to the band index, we shall 
henceforth use them interchangeably. 

In our simulations in figure ll we have chosen g = | — a , which is sufficiently large 
to ensure that the tunnelling probability out of the first band is always appreciable for 
Vo = 1.0 En. With the initial wave function d5|, the condensate starts at t — from the 
state k = (0, 0) in the first zone. For figure [l](a) the acceleration points along the (1, 0) 
direction, therefore according to ^ the condensate would proceed along the k x -axis 
until it meets the zone boundary at k x = -. The wave packet then becomes partly 
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(a) 



(b) 



Figure 2. (a) The first twelve Brillouin zones for the 2D square lattice in the 
extended-zone scheme. The open circles represent reciprocal lattice points, (b) 
shows contours for the corresponding energy landscape when Vo = I.OEr. Here 
k x , k y are both in units of ir/a. For regions not shown here, one can infer the 
zone structure and energy landscape up to the twelfth zone from symmetry of 
the lattice. 



reflected back to the first zone and partly transmitted into the second zone, giving rise 
to the first tunnelling pulse shown in figure flta). The pulse that enters the second zone 
proceeds again according to d9j) until it meets the zone boundary at k = (^,0), where 
three Bragg planes intersect. From figure [21(b) it is clear that the energy gap at this 
point is tiny, thus the pulse can tunnel through entirely. As the wave packet moves 
along, its energy also grows. Eventually the wave packet becomes essentially free, since 
it can overcome any energy gap that it encounters [see (10)]. Notice that this result is 



essentially identical to that for ID optical lattices [TTj [12]. 

For acceleration along the (1,1) direction [see figure [lib)], the fc-space dynamics 
for the condensate becomes more complicated. Again, the condensate starts at t = 
from the state A; = (0, 0) and evolves initially according to (J9l) along the (1,1) direction. 
However, at t — Tb/v2 the state reaches the point k = (-,-), where bands 1 — 4 meet. 
Due to the large energy gap, the condensate can have a partial wave Bragg reflected back 
into the first band [to the point k = (—-,—-)] and a partial wave tunnelling into the 
second band [or the third band, since they are degenerate at this point; see figure[2fb)]. 
For the partial wave entering into the second band, nevertheless, since the acceleration 
g would drive the system to evolve along the (1, 1) direction, the state has to either 
tunnel immediately into the fourth band, or be Bragg reflected within band 2. The four 
atom clouds visible at t = 2.0 Tb in figure [lib) represent each of the above possibilities: 
the atom cloud at the top corner is the partial wave that is Bragg reflected back into 
the first band, executing an oscillatory motion in both real and fc-spaces; the lower one 
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is the partial wave that tunnels directly into the fourth band. The two side pulses are 
partial waves that are Bragg reflected to the points k = (±-,=p-) which subsequently 
evolve under the acceleration g towards the states (— , 0) and (0, — ), respectively. For 
the three partial waves that tunnel out of the first band, the acceleration g continues 
to bring them forward along the (1, 1) direction. The central pulse becomes essentially 
free, as there is no longer any appreciable energy gap along its fc-space trajectory. For 
the two side pulses, however, as one can see from figure [21(b) there remain several gaps 
that they shall cross. 

Since the two side pulses are symmetric, it will suffice to focus our analysis on 
one of them. Let us consider the one on the right. After being Bragg reflected to the 
point (— f , f), the pulse evolves under the pull of g along the (1, 1) direction. It will 
therefore come across an energy gap at k = (^, ||) at the time t = ^^Tb ~ 1.650 Tb. 
For the result in figure fWo), one can rule out any possibility for Bragg reflection at 
this point since the sub-pulses become visible only at t = 2.0 ~ 3.0 Tg (for a more 
detailed check, see later). The next energy gap that the pulse would meet occurs at 
t = -j^Tb — 2.121Tb a t the point k = (^, ^). As we shall explain in greater detail 
below, at this point the pulse turns out to be partially Bragg reflected to the state 
k = (—-, — ) and partially tunnels forward along its original trajectory. Therefore, the 
original pulse is split into two sub-pulses, resulting in the final pattern seen in figure 

lb)- 

To substantiate the arguments in the preceding paragraph, one could consider an 

arbitrary Bloch state on the line segment between the k points (— f , f ) and (^, || ), and 
apply to the state the same acceleration g as above, but now in a controlled manner. 
For instance, starting from the point k = (—f-, f^), we apply the acceleration g for 
the duration At = 1.0 Tb and then turn it off. This would send the pulse to the state 
k = ( ~ 1+ ^ * 5+^2 ^ ^ /oji^ Mi4£) i w hi c h li es between the points f^, ¥■) and 



2 a' 2 a> ~ \ a ' a n 



and then leave it to evolve freely under the lattice potential \2\. We find that in 

ution. If now the 



« 3a ' 3a < 



( — 37T ^ 
\a' a > 

this case the atom cloud remains a single one throughout its time evo 
acceleration is kept on for a longer period of time, so that the pulse could go beyond the 
point k — (-, — ) slightly, one would observe that the two sub-pulses are generated. To 
confirm that Bragg reflection and Landau-Zener tunnelling do occur at (-, — ), one could 
repeat the above procedures but apply now a reduced acceleration along (1,1) direction, 
say, g' = ^^f. One would find in this case only one single pulse appears, even though 
an acceleration duration that could have sent the state beyond k = (-, — ) had been 
applied. What happens here is that the pulse is Bragg reflected from (-, — ) to the 
state (— - , — ), and then tunnels into the next band under the action of g'r\ To further 
endorse our argument, we notice that the directions of the group velocities of the wave 
packets, according to (p]), can be read off from the contours in figure [2|b). Hence the 
sub-pulse that tunnels through from k = (—-,—) has a group velocity almost parallel 

+ There are in fact other possible final states for the Bragg reflection [3 . However, they are suppressed 
here since the optical potential pi) has Fourier components only for Afc = (0,0), (±2n/a, 0), and 
(0,±27r/a). 
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to the (0,1) direction, while the one that tunnels forward from k = (-, — ) would be 
more free-particle like [i.e. its real space trajectory would bend downwards just like a 
projectile; see figure [T|(b)]. From the perspective of energetics, we notice that along the 
(1, 1) direction, tunnelling across the point (— -, — ) would send the wave packet into the 
ninth band (degenerate with the tenth band), while across (-, — ) would bring it into 
the eleventh band (degenerate with the twelfth band). Taking into account the finite 
width of the wave packet, one can find that the former option is evidently energetically 
more favourable. It is also possible to confirm our argument base on Fourier analysis of 
the atom clouds [30] . 

Although a more quantitative analysis of the problem than that presented above is 
possible [T7j, our approach in the foregoing paragraph has its merits. Besides providing a 
check for our argument it also alludes to one important point: in addition to considering 
condensates with zero (center-of-mass) initial velocity (as was done for the simulations 
in figure pj), one could also consider an initial wave function with non-zero wave vector. 
Our analysis thus brings out a useful means for the manipulation of coherent atom 
waves, as we shall now explain. 

3. Manipulation of coherent atom waves 

In order to prepare a condensate with non-zero Bloch wave vector feo, one can start 
from a stationary wave packet with the wave function (J5J). Without turning on the 
optical lattice potential, one applies first an acceleration g (along the k direction) 
to the stationary condensate for the period of time At = -^ e - and then switch it off. 
The condensate would now have acquired the designated wave vector k and the wave 
function would be 



— ex (-? 
2vra 2 6XP l 4<r 



^(P) = \/^T eX P -f^) eX P(^0-p) • (11; 



To project the state onto the lattice eigenstate with Bloch wave vector fco, one could 



switch on the optical lattice and allow the state (11) to evolve for a period of time. 
This would then transfer the condensate into the desired Bloch state. It should be 
noticed that this procedure for state preparation would work better for weak lattice 
potentials (and k a distance away from any energy gaps, of course). For strong lattice 



potentials the Gaussian state ( 11 ) may deviate enormously from the Bloch state k ; thus 
the projection procedure may result in an atom cloud with extremely low density. To 
overcome this difficulty, one may first project the state onto a weak-potential eigenstate 
and afterwards increase the lattice strength steadily. 

In figure [3] we demonstrate manipulation of coherent atom waves utilizing the 
scheme proposed in the previous section. Here we start at t = from the Gaussian 



wave packet (11) with a = 2.5a, N = 2-na 2 , and k = (■£-, ^). The state is left to 



evolve in the optical lattice (|2j) with Vq = 1-0-Er for the period of time At = 1.0 Tg. 
We then turn on the acceleration g = | — a along the (l, l) direction for the duration 
At = 1.0 Tb, sending the state to the point k = (±±§^ f, ^±|^ *) ~ ( i ^, M ^ !L ) i 
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and then switch it off. The condensate is then left to evolve in the lattice potential 
for another period of time At = 2.0 Tg. As one can see in figure [3j the condensate 
begins to split at t = 2.0 Tb and two pulses of atom clouds become visible in the plot for 
t = 3.0Tb- To further split the two pulses, at t = 4.0Tb we apply the acceleration — g, 
which is opposite to that applied earlier. We turn off the acceleration after the duration 
At = 1.0 Tb and allow the pulses to evolve in the lattice freely. The reverted acceleration 
brings the two pulses across the beam-splitting points k — (±-, — ). Therefore, as shown 
in figure |3j at t = 6.0Tb each of the two pulses splits again into two further sub-pulses. 
All these results are consequences of the processes expounded in the previous section. 
One can see that by exploring Bragg reflection and Landau-Zener tunnelling of the 
condensate in accelerated 2D optical lattices, we have furnished a beam splitter that 
may be useful for manipulation of atom waves. By changing the lattice strength Vo 
and/or the magnitude of the acceleration g, one can tune the relative strengths between 
the reflected and the transmitted pulses. Further manipulations of the atom cloud can 
be achieved by extensive exploration of the energy landscape for the optical lattice. For 
instance, one can apply g in arbitrary directions, change g in real time, and/or make 
use of different lattice structures (which may also involve "defects"). A simple example 
for the interference of coherent atom waves is shown in the plot for t = 7.0 Tb in figure 

El 

4. Conclusions and discussions 

The interplay between cold atom physics and solid-state physics has stimulated exciting 
developments in the past decades. In this work, we have studied the dynamics of a cloud 
of ultracold atoms subject to an accelerated 2D optical square lattice by numerically 
solving the GP equation in the regime where mean-field interaction is negligible. We 
have shown that for sufficiently large accelerations the condensate can generate pulses of 
atom clouds that have structures richer than its ID counterparts. Using a semiclassical 
picture, we have demonstrated that these structures can be understood from the k- 
space dynamics of the Bloch state. In particular, Bragg reflection and Landau-Zener 
tunnelling at higher energy bands are shown to be responsible for the pattern generated. 

On the one hand, our result can help envisage 2D atom lasers that have patterned 
output pulses. These patterns are closely connected with the underlying lattice 
structures. On the other hand, we have shown that by appropriately controlling the 
magnitude, direction, and duration of the lattice acceleration, one would be able to 
manipulate the atom wave by exploring the energy landscape of the optical lattice. We 
have demonstrated that a beam splitter for coherent atom waves can be furnished in 
this scheme. 

An important issue that we did not address in this article is the effect of the 
nonlinear mean-field interaction in (fTl). Besides decoherence of the condensate as 
discussed in section 2.1 [121 121], it has been shown that nonlinear effects can cause 
asymmetric Landau-Zener tunnelling and modulation instability for condensates in 
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Figure 3. Numerical simulation for an accelerated condensate that starts from 



the initial wave function (j 1 1 1) with fco = {£~> y~) m the same 2D optical lattice 



as in figure 1 An acceleration g = fxr? is applied along the (1,1) direction in 
a controlled manner as described in the text. As in figure [T] (b), here the time 
scale Tb is given by and the plots are made against the rotated coordinates 
x' and y' in units of the lattice constant a. Notice that for different time frames 
we have shifted the coordinates to different ranges for convenience. 



accelerated 2D optical lattices [201 EI] • We hope to investigate these problems in future 
publications. At the same time, possible applications ensue from this work in the fields 
of atom optics and quantum information sciences are also yet to be explored. 
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